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EVALUATION 


1.  This  is  the  Final  Report  on  the  contract  which  over  the  period  from 

1  December  1976  to  15  March  1980  investigated  the  radiometric  determina¬ 
tion  of  rain  attenuation  at  millimeter  wavelengths,  atmospheric  limitation 
on  the  wideband  millimeter  communication,  backscattering  from  atmospheric 
irregularities,  and  diffusion  of  pulses  in  scatterers.  The  emphasis  was 
placed  on  generating  new  techniques  for  solving  practical  problems. 

2.  The  above  work  is  of  value  since  it  enables  one  to  form  a  quantitative 
understanding  of  the  atmospheric  effect  and  limitation  on  millimeter  wave 
propagation  pertinent  to  the  design  and  operation  of  the  Air  Force  com¬ 
munications  and  detection  systems. 

f\c-  i  t  Ai  ) }  .V/  >  cS 
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Project  Engineer 
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I  Introduction 


In  recent  years,  numerous  studies  have  been  reported  on  the  effects 
of  various  particulate  matter  and  turbulence  in  the  troposphere  on  mille- 
meter  wave  propagation.  In  general,  the  problem  of  wave  propagation  in 
turbulence  has  oeen  well  studied  in  terms  of  the  forward  scatter  approxi¬ 
mation:  This  case  is  characterized  by  the  extremely  small  refractive 
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index  fluctuations  (the  variance  is  of  the  order  of  10  )  and  the  turbu¬ 

lent  eddy  sizes  which  are  much  greater  than  a  wavelength.  In  contrast 
with  the  turbulence  case,  many  practical  situations  present  very  different 
problems.  For  example,  in  the  wave  propagagion  in  fog,  clouds,  or  snow, 
the  refractive  index  of  water  is  very  much  different  from  air  and  the 
particle  sizes  are  comparable  to  or  smaller  than  a  wavelength.  Unfortu¬ 
nately,  there  is  a  definite  lack  of  understanding  of  the  wave  propaga¬ 
tion  and  multiple  scattering  in  the  medium  containing  particulate  matter 
and  other  inhomogenieties  characterized  by  the  particle  sizes  comparable 
to  or  smaller  than  a  wavelength,  particularly  when  the  particulate 
matter  is  densely  distributed.  In  practice,  this  is  the  situation  often 
encountered  for  millimeter  and  optical  waves  in  fog,  clouds,  rain,  hail, 
and  snow.  Also  microwave  scattering  from  rough  terrain  surface  and  the 
volume  scattering  from  terrain  present  the  similar  difficulties  that 
the  inhomogeneities  are  often  of  the  order  of  a  wavelength  and  are 
closely  packed. 

During  the  last  few  years,  we  have  been  concentrating  our  efforts 
on  these  areas  as  indicated  above,  which  have  not  been  clarified  yet. 

We  have  developed  a  diffusion  theory  which  may  be  applicable  to  dense 
random  media.  However  the  range  of  validity  of  this  approximation  has 
not  been  established  yet.  At  present,  therefore  we  have  not  obtained 
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complete  formulas  for  angular  (field-of-view)  broadening,  pulse  broaden¬ 
ing,  amplitude  and  phase  fluctuations,  temporal  and  spatial  spectra, 
polarization  effects,  the  effects  of  beam  waves,  the  effects  of  size 
distributions  ana  the  densities  when  millimeter  waves  are  incident  on 
various  media  including  the  particulate  matter  in  the  atmosphere  and 
snow,  vegetation  and  other  random  media  in  the  terrain.  We  hope  to 
continue  to  develope  some  fundamental  formulations  for  these  problems 
and  to  obtain  useful  numerical  predictions  for  the  wave  propagation  and 
scattering  in  these  media. 

1 1  Radiometric  Determi nation  of  Rain  Attenuation  of  Millimeter  Wavelength 

In  rain  attenuation  measurement  using  a  radiometer,  scattering  is 
usually  negligible  compared  with  absorption  under  30  GHz  and  light  rain. 

At  mi  1 1 imeter  wavelength ,  however,  it  is  necessary  to  include  multiple 
scattering  effects.  We  have  conducted  a  study  of  the  effects  of  scattering 
taking  into  account  the  polarization  effect  and  the  drop  size  distribution. 

The  formulation  is  based  on  the  equation  of  transfer  with  Stokes 
parameters.  We  assumed  that  rain  is  located  in  a  plane-parallel  region 
and  that  the  thermal  radiation  from  rain  droplets  obeys  the  Kirchhoff  law. 
The  absorption  and  scattering  cross  sections  and  the  phase  functions  are 
calculated  using  the  Mie  solution  and  taking  into  account  the  Laws  and 
Parsons  size  distribution.  The  index  of  refraction  of  water  is  calculated 
by  using  the  formula  given  by  Saxton. 

For  a  plane-parallel  problem,  the  specific  intensities  are  functions 
of  height  and  angle  only.  Letting  Ig  and  Ir  be  the  vertical  and  horizontal 
(parellel)  components  of  the  specific  intensity,  the  equation  of  transfer 
becomes 
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All’  Aj2*  ^21*  ^22  are  9’ven  ^  Zekera.  The  general  solution  is  a  sum 
of  a  particular  solution  and  a  complementary  solution.  The  complementary 
solution  can  be  obtained  by  the  eigenvalue  technique. 

The  general  boundary  condition  at  ground  is  given  by 

T(s)  =  [1  -  /  R(s,i)p.ds^Tg  +  /  R(s,i  )T(i  )y.dsi. 

where  s  and  i  are  unit  vectors  in  the  directions  of  scattering  and 
incidence,  and  the  integrals  are  taken  over  the  upper  hemisphere.  If  the 
diffuse  scattering  is  uniform,  then  R(s,i)  is  constant  and  is  equal  to 
(Albedo)/n,  yielding 

T(s)  .  [1  -  Rg]Tg  *  2Rg  /*  Kiliylu, 


where  R  is  the  albedo  of  the  surface. 

g 

The  general  equation  of  transfer  was  solved  using  the  eigenvalue 
technique.  Calculations  are  made  for  the  frequency  ranges  from  30GHz  to 
120GHz  and  various  precipitation  rates.  For  example,  at  30GHz,  for  the 
rain  temperature  of  273°K,  the  precipitation  rate  of  12.5nim/hr,  the  ground 
temperature  of  283°K,  and  the  rain  layer  thickness  of  2km,  the  calculated 
temperature  for  vertical  and  horizontal  components  are  found  to  range  from 
171°K  to  260°K  and  171°K  to  257°K  respectively  as  the  observed  direction 
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varies  from  vertical  to  horizontal.  But  at  1206Hz,  the  observed  temperatures 
for  vertical  and  horizontal  polarizations  vary  from  248° K  to  266°K  and  248°K 
to  263°K  respectively.  Based  on  these  calculations,  the  differences  between 
the  true  rain  attenuation  and  the  attenuation  calculated  from  the  temperature 
measurements  and  the  assumed  rain  temperatures  can  be  calculated. 

The  computer  calculations  are  made  fro  30  GHz,  60  GHz,  and  120  uHz  and 
the  rain  layer  thicknesses  of  1  km,  2  km,  and  3  km.  The  albedo  of  ground 
is  also  varied,  but  the  effect  seems  to  be  small.  The  precipitation  rates 
used  include  0.25  mm/hr,  12.5  mm/hr,  and  100  mm/hr. 

With  the  radiometric  measurement  technique,  the  attenuation  by  rain 
is  usually  estimated  by  the  apparent  attenuation 

T  (1) 

A  (db)  =  10  log  T-S— 
m  t,r 

where  T  is  the  rain  temperature  and  T  are  the  vertical  and  horizontal 
brightness  temperatures ,  respecti vely .  This  approximation  is  only  true  for 
frequencies  below  30  GHz.  At  frequencies  above  30  GHz  and  heavier  pre¬ 
cipitations,  Eq.  (1)  is  no  longer  valid  because  the  scattering  effects 
become  comparable  with  absorption.  The  true  attenuation  should  be  calculated 
in  the  following  way: 

po  .z 

At(db)  -  (10  log  e)^-  (2) 

where  zq  is  the  thickness  of  the  rain  layer  and  e  is  the  observation  angle. 

The  true  attenuation  A^.  is  found  to  be  greater  than  the  apparent  attentuations 
(A.  A  )  and  their  differences  have  been  calculated.  For  a  small  pre- 

x.  f 

cipitation  of  0.25  mm/hr,  the  differences  between  A^  and  A^,  Ar  are  under 
3  db  for  frequencies  below  120  GHz  and  Eq.  (1)  gives  a  good  approximation. 
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The  differences  increase  as  the  frequency  and  precipitation  increase 
because  more  scattering  occurs  with  higher  frequencies  and  heavier 
precipi tations . 

It  is  also  useful  to  plot  the  true  attenuation  (A$)  vs.  the  vertical 
apparent  attentuation  (A^)  so  that  the  true  attenuation  can  be  conveniently 
found  from  the  known  vertical  apparent  attenuation.  We  compared  our  results 
with  Zavody's  theoretical  study  and  Chadha-Lane ' s  experimental  data. 

Complete  descriptions  of  our  studies  are  included  in  the  M.S.  thesis  by 
R.  Cheung.  A  paper  summarizing  this  study  has  been  published.  Our  study 
and  Zavody’s  are  both  based  on  spherical  rain  droplets.  Chadha-Lane's 
experiment,  however,  showed  noticeable  differences  between  horizontal  and 
vertical  polarizations,  indicating  the  effect  of  non-spherical  rain  droplets. 

Ill  Atmospheric  Limitation  on  Wide  Band  Communication 

We  have  conducted  a  study  of  transmission  characteristics  of  millimeter 
and  optical  waves  under  various  atmospheric  conditions.  We  have  concentrated 
our  attention  to  the  atmospheric  turbulence,  rain  and  fog.  We  can  summarize 
the  general  characteristics  as  follows: 

(a)  Millimeter  Wave 

Turbulence:  The  scale  sizes  of  the  turbulence  are  so  much  greater  than  a 
wavelength,  and  the  angular  spread  of  the  wave  is  very  small.  Also  the 
absorption  is  much  smaller  than  the  scattering.  Therefore,  the  forward 
scatter  theory  is  applicable  in  almost  all  ranges  of  practical  interest 
and  the  theory  is  well  developed. 

Rain:  The  droplet  sizes  are  comparable  to  a  wavelength  and  the  scattering 
cross  section  is  of  the  same  order  as  the  absorption  cross  section.  We 
have  conducted  extensive  studies  on  this  problem.  We  have  examined  the 
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range  of  validity  of  the  first  order  multiple  scattering  solution  in 
terms  of  the  precipitation  rate  and  the  propagation  distance.  Above 
10  GHz,  the  multiple  scattering  due  to  rain  may  become  significant 
depending  upon  the  precipitation  rate,  the  thickness  of  the  rain  layer, 
and  the  observation  angle.  We  considered  the  problem  of  a  wave  normally 
incident  upon  a  plane-parallel  rain  region.  The  scattered  wave  is 
observed  on  the  ground  as  a  function  of  the  elevation  angle  and  polar¬ 
izations,  and  it  is  calculated  for  10,  30  and  120  GHz  at  the  precipitation 
rate  of  12.5  mm/hr  with  the  rain  layer  thicknesses  of  1  km  and  3  km. 

The  scattering  characteristics  of  raindrops  are  calculated  using  the  Mie 
solution  and  the  Laws-Parsons distribution.  Both  the  first  order  multiple 
scattering  theory  and  the  radiative  transfer  theory  are  used  and  the 
results  are  compared.  It  is  shown  that  at  10  GHz,  these  two  methods  yield 
almost  identical  results.  However  at  30  GHz  and  120  GHz,  the  first  order 
scattering  is  smaller  than  that  obtained  from  the  radiative  transfer 
theory,  indicating  the  effects  of  higher-order  scattering.  It  is  also 
shown  that  the  horizontal  polarization  is  in  general  greater  than  the 
vertical  polarization  due  to  the  scattering  characteristics  of  raindrops. 

As  the  angle  from  the  zenith  increases,  the  scattering  in  vertical  polar¬ 
ization  decreases  at  10  and  30  GHz.  However  at  120  GHz,  the  scattering 
in  both  polarizations  decreases  with  the  angle,  indicating  the  effect  of 
forward  scattering.  Based  on  the  above  study,  the  incoherent  (fluctuating) 
intensity  for  a  receiver  with  a  given  field  of  view  or  the  receiving  angle 
can  be  calculated.  For  example,  with  the  receiver  beamwidth  of  5°  and 
the  rain  thickness  of  3  km,  the  incoherent  intensity  is  approximately  50  db 
down  from  the  coherent  intensity  at  10  GHz,  but  it  is  20  db  down  at  120  HGz. 
The  incoherent  intensity  for  the  rain  thickness  of  1  km  is  about  5  db  less 
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than  that  for  the  rain  thickness  of  3  km. 

Fog :  The  fog  particle  sizes  are  of  the  order  of  several  microns  and  the 
Rayleigh  scattering  is  applicable.  The  first  order  multiple  scattering 
theory  should  be  applicable  in  this  case. 

(b)  Optical  Wave 

Turbulence:  The  scale  sizes  are  much  greater  than  a  wavelength  and  the 
scattering  is  much  greater  than  the  absorption.  The  well  developed  forward 
scatter  theory  is  applicable. 

Rain:  The  droplet  sizes  are  still  greater  than  a  wavelength  and  the  scatter¬ 
ing  is  greater  than  the  absorption.  The  forward  scatter  theory  should  be 
appl icable  here. 

Fog:  The  particle  sizes  are  comparable  to  a  wavelength  and  the  scattering 
is  much  greater  than  the  absorption. In  this  case,  the  diffusion  should  be 
predominant.  We  have  conducted  extensive  studies  on  the  transmission 
characteristics.  We  compared  the  exact  solution  with  the  diffusion  solution 
to  determine  the  range  of  validity  of  the  diffusion  theory. 

(c)  Coherence  bandwidth  of  scattering  medium 

The  coherence  bandwidth  of  the  atmospheric  turbulence  and  hydrometeors 
is  an  important  parameter  in  communication.  We  have  made  an  extensive  study 
on  two-frequency  mutual  coherence  function  and  pulse  characteristics  in  fog, 
cloud,  and  turbulence.  For  particles,  whose  sizes  are  large  compared  with 
a  wavelength,  we  obtained  the  coherence  bandwidth  wco^  as  follows: 
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-  — “ -  when  p  o  z  >>  1 

2  ns 

L>  o  z- 

K  n  s 

where Pn  is  the  number  density,  o$  is  the  scattering  cross  section,  and  the 
scattering  amplitude  f  is  approximated  by 

M2  =  °s  exp(-ups2) 

The  pulse  shape  for  a  delta  function  input  for  p^z  >>  1  is  given  by 

TT  O)  ,  00  x 

G(t)  =  [_coh]  l  (2n  +  l)(-l)n  exp{-  [(2n  +  l)Jf2  T} 

4  n=o 


v/here  T  =  ^oh( t  -  f)  . 

For  turbulence,  we  obtained  the  coherence  bandwidth: 

0.31  u 

_  _ O 

“coh  "  12/5 

a 

\ 

2 

where  w  is  the  carrier  frequency  and  a  is  the  log-amplitude  variance  for 
°  X 

the  Rytov  solution.  We  also  obtained  a  universal  pulse  shape  in  the  strong 
fluctuation  region. 

(d)  Forward  Scatter  Theory  and  Diffusion  Theory 

The  scale  sizes  of  atmospheric  turbulence  are  usually  so  large  compared 
with  a  wavelength  that  the  angular  broadening  is  small  and  the  forward-scatter 
theory  is  applicable.  For  waves  in  scatterers  such  as  clouds  and  fog,  the 
particle  sizes  are  comparable  to  a  wavelength  and  the  diffusion  phenomena 
become  dominant  at  relatively  short  optical  distances.  It  is  therefore  • 

necessary  to  examine  the  two  distinct  states  of  propagation  and  scattering 
in  random  media:  forward  scatter  and  diffusion.  In  forward  scatter  region 
backscattering  is  negligible  compared  with  forward  scattering  while  there 
is  as  much  backscattering  as  forward  scattering  in  diffusion  region.  The 
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forward  scatter  theory  is  based  on  parabolic  approximation  where  the 
scattering  is  confined  in  a  small  forward  angular  region.  In  contrast, 
the  diffusion  theory  is  based  on  an  approximation  where  the  angular 
spectrum  is  almost  uniform.  They  are  two  asymptotic  theories  which  can 
be  derived  from  a  general  radiative  transfer  equation.  Some  formulations 
and  approximation  solutions  in  these  two  regions  for  both  cw  and  pulse 
cases  have  been  obtained.  For  a  1 ine-of-sight  problem,  the  diffusion 
solution  is  seen  to  be  applicable  for  optical  scattering  distance  of 
greater  than  one.  The  transmitted  flux,  the  angular  spectrum  and  the 
pulse  broadening  are  obtained.  For  backscatteri ng,  the  diffusion  solution 
is  still  applicable  except  for  a  short  initial  time.  This  is  to  be  expected 
as  the  diffusion  theory  is  not  applicable  near  the  boundary.  The  relation¬ 
ships  between  the  forward  scatter  theory  and  the  diffusion  theory  are  not 
completely  clarified  yet  and  their  range  of  validity  should  be  established. 


I V  Effects  of  Transmitter  and  Receiver  Characteristics  on  the  Wave 
FI uctuation. 

We  have  made  extensive  study  on  the  effects  of  transmitter  and  receiver 

characteristics  on  amplitude  and  phase  fluctuations  in  the  weak  fluctuation 

2  2 

region.  The  variances  of  amplitude  <x  >  and  phase  <S  >  fluctuation,  including 
transmitter  and  receiver  characteristics,  are  given  by 

po_L  1 


•x2>  = 


^/d„ 


<S2>  = 


po  L  1  ,  , 

t5-  / +  Re  rr2 


Q  .".td.y.V  ! 

1  u  '  "  a  n  C 


where 
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a  ,  pa  L-i2kL  > 

Q,  =  -±  (1-n)2  +  1  +  — ! l- - 

c  a_  a_  4 


where 


.  _  2.77  _  _  2.77  _  _  2.77 


and  etb,  @rb  and  e  b  are  the  half-power  beamwidths  of  the  transmitter, 

2 

receiver  and  scatterer,  respecti vely.  We  have  calculated  both  <x  > 

2 

and  <S  >  for  a  rain  medium  with  precipitations  of  12.5  mm/hr  and  50  mm/hr 
at  30,  60,  90  and  120  GHz  and  at  optical  frequency  x  =  0.6  y. 

At  microwave  frequencies  (30,  60,  90  and  120  GHz),  <x2>  and  <S^>  are 
very  close  to  each  other.  They  both  increase  with  the  distance  until  they 
saturate  at  large  distance.  They  also  increase  with  the  increase  of 
frequency  and  precipitation.  At  optical  frequency  (  X  =  0.6  y,  p  -  12.5  mm/hr), 
the  difference  between  <x2>  and  <S2>  can  be  significant. 

These  studies  are  made  using  the  Gaussian  model  for  the  scattering 
pattern.  However,  further  studies  are  needed  using  more  realistic 
scattering  patterns.  It  is  also  necessary  to  extend  the  results  to  the 
strong  fluctuation  region. 


V  Backscattering  of  a  Pulse 

We  have  conducted  a  study  of  backscattered  pulse  in  two  asymptotic 
cases:  forward  scatter  and  diffusion. 

The  cummulative  forward-scatter  single-backscatter  (CFSB)  approxi¬ 
mation  were  discussed  by  deWolf  (IEEE  Trans.,  AP-19,  254-262,  1971)  and 
this  is  applied  to  a  plane  wave  incident  upon  a  slab  for  calculation  of 
the  backscattered  two-frequency  mutual  coherence  function,  from  which  the 
backscattered  pulse  shape  is  found  by  Fourier  transformation.  The  response 
to  a  general  pulse  may  then  be  obtained  by  convolution.  To  employ  this 


v* 
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technique,  the  forward-propagating  two- frequency  mutual  coherence  function 

is  expressed  as  the  solution  to  a  parabolic  equation.  Multiple  scattering 

effects  are  included,  under  the  restrictions  of  small  angular-spread  of  the 

wave  and  a  narrow- bandwidth  envelope  for  the  incident  pulse.  The  governing 

parameters  for  the  pulse  shape  are  the  depth  of  the  medium  measured  by  the 

optical  distance  x,  the  albedo  Wq  and  a  parameter  which  is  inversely 

related  to  the  angular  width  of  the  averaged  scattering  pattern  of  the 

medium.  When  the  scattering  pattern  is  a  Gaussian  function,  the  latter 
2  2 

parameter  is  y  =  i  kQ  /4,  where  a  is  the  scale  length  of  the  correlation 

function  of  refractive  index  fluctuation  and  k  is  the  wave  number.  The 

o 

range  of  validity  of  the  CFSB  approximation  is  found  to  be  approximately 
2 

t  s  (26q/3)  >/Wq  where  0Q  is  the  allowable  angular-spread  within  the  medium. 

Although  the  overall  temporal  pulse  width  is  always  approximately  2L/c, 

where  L  is  the  actual  depth  of  the  slab  and  c  is  the  velocity  of  propagation, 

the  pulse  shape  is  best  preserved  when  the  time  scale  is  normalized  to 

2L/ct.  The  effect  of  non-unity  albedo  is  to  hasten  attenuation.  The 

effect  of  y  is  surprisingly  weak,  as  drastic  increases  in  y  serve  only 

to  sharpen  the  leading  and  trailing  edges  of  the  overall  pulse  shape. 

Thus,  asymptotic  solutions  may  be  valuable,  and  a  closed-form  solutions 

as  y  -*■  °°  is  presented.  When  scattering  is  instead  characterized  by  a 

Kolmogorov  process,  the  characteristic  parameter  k  LQ  is  introduced,  where 

Lq  is  the  outer  scale  length.  The  resulting  pulse  shape  differs  greatly  from 

the  pulse  shape  for  a  Gaussian  spectrum.  I’  this  case  the  range  of 

validity  for  the  CFSB  approximation  is  found  to  be  x  <  (kQLo)^30o^/12Wo. 

The  effect  of  changing  k  L  is  quite  pronounced;  however  a  closed-form 

o  o 

asymptotic  solution  for  the  pulse  shape  as  kQL0  is  found,  which  in 
extreme  conditions  of  short  wavelength  of  a  large  otter  scale  size  may  be 


••  >  ^ 
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a  good  approximation  to  actual  behavior. 

Diffusion  of  a  short  pulse  in  a  dense  scattering  random  medium  has 
previously  been  studied  using  the  time-dependent  equation  of  transfer 
and  its  diffusion  approximation  (A.  Ishimaru,  JOSA,  68,  1045-1050,  1978). 
Ue  studied  a  diffusion  solution  for  a  beam  wave  incident  upon  a  slab  of 
a  dense  scattering  medium.  The  backscattered  pulse  shapes  are  calculated 
from  the  theory  and  compared  with  the  experimental  data.  The  experimental 
system  was  constiucted  by  modifying  the  laser  range-gating  system 
(A.  P.  Bruckner,  Appl.  Opt.,  J7,  3177-3183,  1978).  Using  high  power  laser 
pulses  (^500  MW/cm  )  and  an  ultrafast  shutter  (CS^  between  two  crossed 
polarizers),  we  were  able  to  record  backscattered  pulse  shapes  in  a 
picosecond  ranye.  A  beam  of  pulses  with  the  pulse  width  of  approximately 
10  picoseconds  at  A  =  0.53  pm  is  incident  upon  an  aqueous  solution 
containing  latex  microspheres.  The  latex  spheres  have  diameters  ranging 
from  0.5  to  50  pm  and  their  concentrations  are  in  the  range  of  0.1  to  30% 
in  weight.  The  backscattered  pulse  shapes  are  digitized  and  stored  for 
statistical  processing.  The  pulse  broadening  for  a  10%  solution  is  a 
few  picoseconds,  while  the  broadening  for  a  1%  solution  is  approximately 
5  to  10  picoseconds.  The  pulse  broadening  increases  as  the  concentration 
is  decreased  and  for  a  0.1%  solution,  the  pulse  broadening  extends  to 
several  tens  of  picoseconds.  The  general  pulse  shape  is  characterized 
by  a  sharp  rise  time  and  slow  tail.  When  the  scatterers  are  absorbing, 
the  tail  part  tends  to  suffer  attenuation.  The  diffusion  solutions  are 
calculated  using  parameters  obtained  from  the  Mie  solution  and  the  known 
concentration.  It  is  demonstrated  that  the  theoretical  predictions  and 
the  experiments  agree  well  not  only  in  the  pulse  shape  and  the  broadening, 
but  in  the  relative  magnitudes  of  the  pulse  height  for  different  particle 
sizes  and  concentrations  as  well. 


VI  Diffusion  of  a  Beam  Wave  in  a  Scattering  Medium 


The  diffusion  solution  of  a  collimated  beam  wave  normally  incident 
upon  a  plane  parallel  layer  of  scatterers  have  been  examined  making  use 
of  the  diffusion  equation  and  boundary  conditions  given  in  Chapter  9  of 
Ishimaru's  "Wave  Propagation  and  Scattering  in  Random  Media"  -  Vol .  1. 
By  applying  Green's  function  technique  and  using  the  Hankel  transform, 


the  solution  for  the  Average  Diffused  Intensity  has  been  found: 

w02  i  c  c 

u  (r,z)  =  f  ^  J  (Ar)'*W  4  X  T  A(z)  +  4FFT  B(z)  ’ 


where  f  =  (a2+K^?)2,  is  the  beam  size  of  the  incident  collimated  wave 
with  Gaussian  profile,  and  A(z)  and  B(z)  are  known  functions  of  z. 


The  beam  spread  W 


P 


is  defined  as: 


2ii  /  F  (r,d)  r  d  r 
o  z 


where  F^  is  the  normal  component  of  the  transmitted  flux  at  the  exit  end 
of  the  layer. 

For  example,  at  A  =  5pm,albedo  w  =  0.8758,  m  =  0.7774,  WQ  =  1  cm, 

W  - 

the  quantity  (J-)  increases  to  10'5  when  the  optical  thickness  t  of  the 
O  5 

layer  is  10,  and  becomes  7.0x10  when  t  =  20.  We  intend  to  continue  to 
study  the  transmitted  and  reflected  flux  of  the  layer  and  the  diffused 
flux  distribution  F.z  (r,d)  at  the  end  of  the  layer.  The  results  for  the 
small  optical  thickness  are  compared  with  the  forward  scattering  results. 
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